Protection of quantum systems by nested dynamical decoupling 
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Based on a theorem we establish on dynamical decoupling of time-dependent systems, we present a scheme 
of nested Uhrig dynamical decoupling (NUDD) to protect multi-qubit systems in generic quantum baths to ar- 
bitrary decoupling orders, using only single-qubit operations. The number of control pulses in NUDD increases 
polynomially with the decoupling order. For general multi-level systems, this scheme can preserve a set of 
unitary Hermitian system operators which mutually either commute or anti-commute, and hence all operators 
in the Lie algebra generated from this set of operators, generating an effective symmetry group for the system 
up to a given order of precision. NUDD can be implemented with pulses of finite amplitude, up to an error in 
the second order of the pulse durations. 
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I. INTRODUCTION 



Both in high-precision magnetic resonance spectroscopy 
and in quantum computing, it is essential to suppress un- 
wanted couplings within a quantum system and between the 
system and its environment (or bath). Such couplings result in 
population relaxation, phase randomization (pure dephasing), 
and more generally, unwanted evolution of certain system 
operators. Dynamical decoupling (DD) is a so-called open- 
loop control scheme to average out the system-bath interac- 
tions through stroboscopic operations of the system (with- 
out direct control on the environment). DD originated from 
the Hahn echo lit] and has evolved into variations with com- 
plicated sequences |01 for high-precision d, H] and multi- 
dimensional [Si] magnetic resonance spectroscopy. When the 
field of quantum computing was opened up, DD was intro- 
duced to protect qubit coherence M 01- By unitary sym- 
metrization procedure Ull^, DD cancels errors of quantum 
evolutions up to the first order in the Magnus expansion, and 
the corresponding cyclic scheme is referred to as periodic DD. 
A geometric understanding of the symmetrization procedure 
was given in Ref . ifioll . To eliminate errors to the second or- 
der in the Magnus expansion, mirror-symmetric arrangement 
of two DD sequences (SDD) can be used |@]. A particularly 
interesting scheme is the concatenated DD (CDD) ifTll - fTill . 
which uses recursively constructed pulse sequences to elim- 
inate decoherence to an arbitrary decoupling order (defined 
as the power of the total evolution time, which is assumed 
short). The performance of CDD was experimentally demon- 
strated for spins in solid-state environments (Tsl). The number 
of pulses used in CDD, however, increases exponentially with 
the decoupling order. Since errors are inevitably introduced in 
each control pulse in experiments, finding DD schemes with 
fewer control pulses is desirable. 

For suppressing pure dephasing of single qubits (two-level 
systems) subjected to unidirectional noises, a remarkable ad- 
vance is the optimal DD discovered by Uhrig lfl6ll in a spin- 
boson model. Uhrig DD (UDD) is optimal in the sense that the 
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number of control pulses is minimum for a given decoupling 
order. It was later conjectured ITtI [TsIl and then rigorously 
proven lfl9ll that UDD is model-independent for any two-level 
systems coupled to a finite quantum bath. It was shown that 
UDD also works for suppressing longitudinal relaxation lfl9ll . 
The ideal (5-pulses assumed in DD can be generalized to in- 
clude some components of finite amplitude lfl9ll . Recently, a 
method to incorporate shaped pulses of finite amplitude into 
UDD paves the way of realistic experiments ||20tl. UDD was 
first verified in experiments by microwave control of trapped 
ions in various artificial classical noises ll2ll - t23ll , and then 
UDD against realistic quantum noises was realized for radi- 
cal electron spins in irradiated malonic acid crystals ll24ll . 

For suppressing the general decoherence of single qubits 
(including both pure dephasing and longitudinal relaxation), 
the concatenation of UDD sequences (CUDD) was proposed 
to reduce the number of control pulses jzsll . For suppressing 
the decoherence up to an order A^, the number of pulses re- 
quired in CUDD is ~ (A^ + 1)2'^, which is considerably less 
than ~ 4'^ as in CDD. Recently, West et al proposed a much 
more efficient scheme, called quadratic DD (QDD), to com- 
bat general decoherence of a qubit |2^. QDD is constructed 
by nesting two levels of UDD sequences, using (A^ -i- 1)^ con- 
trol intervals to achieve the A^th decoupling order Numerical 
search indicates that QDD is near-optimal as it differs from the 
optimal solutions by no more than two pulses for a small de- 
coupling order (A^ < 4) ll26ll . The validity of UDD can be 
extended to analytically time-dependent Hamiltonians ll27ll . 
This extension seems to validate QDD since the UDD se- 
quence on the outer level can be viewed as acting on a time- 
dependent Hamiltonian resulting from the UDD control on the 
inner level. However, as we will show in Sec. |III1 the effec- 
tive Hamiltonian resulting from the inner UDD sequences is 
only piecewise analytic in time, and actually there are counter 
examples in which UDD on the outer level does not achieve 
the designed decoupling order if the order of UDD on the in- 
ner level is odd and lower than the order of the outer UDD. 
Thus it remains an open question why QDD works. In the at- 
tempt to prove the validity of QDD, we establish a theorem: 
UDD applies to time-dependent Hamiltonians, regardless of 
their analytic properties, as long as the Hamiltonians between 
two adjacent pulses of the UDD sequence are symmetric and 
have the same function form of relative time between the ad- 
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jacent pulses (see Sec. Ill Pi . Therefore we give a proof of the 
vaUdity of QDD with even order UDD on the inner level. We 
note that the validation of QDD is still incomplete since the 
theorem mentioned above does not apply to QDD with odd 
order UDD on the inner level. 

So far, the research of optimal DD focuses on the single- 
qubit decoherence problem, with some attempts on optimal 
DD to multi-level systems with prior knowledge of the ini- 
tial states ||28[ l29ll . For practical large-scale quantum com- 
puting, the issue is the decoherence, or more generally, the 
decay of quantum correlations (such as entanglement) of cou- 
pled multi-qubit systems. Therefore it is highly desirable to 
have a general arbitrary-order DD scheme for multi-qubit sys- 
tems, with the number of pulses as small as possible. In addi- 
tion, the DD scheme should preferably involve relatively easy 
implementations, for example, single-qubit operations. Re- 
alizations of the symmetrization procedure lH |2I by averag- 
ing over so-called nice error bases ifsoll are explicitly given 
in Ref. llsill . However, each control operation in general is 
comphcated and may involve manipulation on all qubits for 
multi-qubit systems. It would be of practical interest if there 
is an explicit, systematic, and efficient way to protect a par- 
ticular set of operators of a quantum system to an arbitrary 
order 

In this paper, we give systematic and expUcit DD schemes 
to protect multi-qubit systems arbitrarily coupled to quantum 
baths. The schemes are realized by protecting a set of mutu- 
ally commuting or anti-commuting unitary Hermitian system 
operators on different levels. We call it a mutually orthog- 
onal operation set (MOOS) in this paper For example, the 
Pauli operators of qubits form an MOOS. The inner levels 
of DD control of operators in an MOOS are not affected by 
the outer levels of DD control. Furthermore, based on the 
above-mentioned theorem on DD of time-dependent systems, 
higher order protection of an MOOS can be achieved by nest- 
ing even-order UDD sequences on different levels. If a set of 
system operators is protected by such nested UDD (NUDD), 
then all system operators in the Lie algebra generated from 
this set of operators are protected to the same decoupling or- 
der, which indeed generates an effective symmetry group iH] 
of the system up to an error of the decoupling order. For 
multi-qubit systems, each control operation in DD only in- 
volves single-qubit manipulation. In addition, we will show 
that NUDD can be implemented with pulses of finite ampli- 
tude, which approximate ideal 5-pulses up to an error in the 
second order of the pulse durations, with the same pulse shap- 
ing as in Ref. ll32ll . For a general multi-level quantum system, 
we can also construct an MOOS and use NUDD to generate an 
effective symmetry group to a given decoupling order It can 
be shown, however, that for a general M-level system, there 
may exist no MOOS to generate the whole basis of su(M) 
algebra and hence the whole SU(M) symmetry group. Fur- 
ther research is still needed to design efficient DD schemes 
(as compared with CDD) to protect general multi-level sys- 
tems to higher orders. 

This paper is organized as follows. In Sec. [Ill we present 
a general theory on protection of an MOOS by DD; NUDD 
is given based on a theorem established for UDD on time- 



dependent systems. In this Section, the pulses in DD are 
assumed instantaneous. In Sec. |III] we discuss NUDD on 
multi-qubit systems. In Sec. |IV] we discuss DD with finite- 
amplitude pulses. Finally, we draw the conclusions in Sec. |V] 

II. PROTECTION OF SYSTEM OPERATORS 

A. General formalism: MOOS 

We consider a quantum system coupled to a general finite 
quantum bath, with a time-independent Hamiltonian 

H^Hs+Hb + Hsb, (1) 

where Hs is the system Hamiltonian, Hb the bath one, and 
Hs B the system-bath interaction. We aim to find a sequence 
of stroboscopic operations 0\, O2, ■ ■ •, Og at times Ti, To, ■ ■ ., 
Tq in increasing order so that after the controlled evolution 
from f = to f = r, a set of system operators 

are conserved up to an error of O (r'*'^'), i.e., for {% - 1) 

U = ■ ■ ■ 026-'"^''^--^^ ^Oie-'"^^ , (2) 

we have 

f/^e;C/ = e/ + o(r^+'). (3) 

In this paper we assume that the bath is bounded in spectrum 
so that the perturbation expansion of [/ in a short time T is 
possible. To make DD efficient, it is required that the evolu- 
tion of the system induced by the decoupling field be faster 
than the unwanted dynamics |@]. Here we assume that the 
control pulses are instantaneous and arbitrarily strong. Using 
pulses of finite duration and finite amplitude will be discussed 
in SecHVl Preferably, for a multi-qubit system, the operations 
Oj should contain only single-qubit operations. We also wish 
to use as few as possible operations to achieve a given order 
(A^) of decoupling precision. 

We note that when a set of operators [Qj] is preserved to a 
certain decoupling order, then all operators obtained by com- 
mutation i[Qj,Qk], anti-commutation [Qj,Qk]+ = QjQk + 
QkQj, linear combinations, and their repetitions are also pro- 
tected to the same decoupling order All these protected oper- 
ators form a Lie algebra, since they form a linear vector space 
and are close under the commutation operation. The Lie alge- 
bra defines a dynamically generated effective symmetry group 
Sq of the system up to the decoupling order 

In particular, let us consider the protection of a set of oper- 
ators {Qy) in which each pair of elements either commutes or 
anti-commutes. In our schemes, it is required that the opera- 
tors {Qj) be unitary and satisfy = +L Thus we choose Qy 
to be unitary and Hermitian, i.e., 

q2 = Q"';q, = L (4) 

Note that is a parity kick operator described in Ref. Issll . 
We will use the unitary Hermitian property of Hj to construct 
DD for protection of Q^. 
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Definition 1. An MOOS is defined as a set of operators which 
are unitary and Hennitian and have the property that each 
pair of elements either commutes or anti-commutes. 

The commutation property of operators in an MOOS is im- 
portant for constructing higher-order DD schemes, via the fol- 
lowing theorem. 

Tlieorem 1. For two operators Q.i and Q2 either commuting 
or anti-commuting with each other, if a unitary evolution U{T) 
during a short time T commutes with Qi up to an error of 
0(7^^+1), then ^2U{T)Q.2 also commutes with D.\ up to an 
error of 0{t^^^). 

Proof. We directly calculate the commutator 

[Qi,Q2t/(r)Q2] = (-i)''Q2 [Qi,f/(r)]Q2 = 

where 77 - +1 (-1) for Qi commuting (anti-commuting) with 
^2- □ 

With this theorem, certain DD sequences protecting a given 
operator in an MOOS can be used as units to construct an 
outer level of DD protection of another operator in the MOOS, 
without affecting the DD effect of the inner level control. For 
operators which do not form an MOOS, the outer level control 
in general may interfere with the inner level control. For ex- 
ample, let us consider Q.i = cr^ and Q2 = (ctv + cr,,^ / V2 (the 
Pauli matrix along the direction equally dividing the angle be- 
tween the X- and y- axes). Suppose a DD sequence protects 
(Tv as U{T) = g-ia-xT+0(T ) jjjgjj apply a Hahn echo to 
protect O2, by the following evolution 

U2(2T) = ^2U(T)n2U(T) = e--,r+o(r-)^-,..,T+o(r-)_ 

After the control on the outer level, the evolution actually does 
not commute with cT v even in the leading order of T. 

As we will show in Sec. IIIBI applying an element of 
an MOOS corresponds to a symmetrization procedure IH |^ 
on one level, and protecting an MOOS will iteratively sym- 
metrize the system evolution. An MOOS itself does not form 
a group or an algebra. However, it can generate a certain Lie 
algebra by commutation, anti-commutation, linear combina- 
tions, and repetitions. Therefore the protection of all operators 
in this Lie algebra is realized by the protection of this MOOS. 
For example, in a single-qubit system, pulses of 7r-rotation 
control of the qubit about the z-axis, when arranged according 
to the UDD timing, protect the Pauli matrix cr- and dynam- 
ically generate the symmetry group SO(2) or U(l), up to an 
error of 0(T^^^) |[l9ll . For a pure dephasing Hamiltonian, 
which has U(l) as its intrinsic symmetry group, the system 
under the UDD protection has the full SU(2) symmetry. An- 
other example is the Mh order QDD sequence ll26ll consisting 
of pulses of cTj and cr,,, which protects cr^ and cr,, and hence 
all operators in the su(2) algebra (including cr~), and therefore 
the qubit has a dynamically generated symmetry group SU(2), 
up to an error of O (t^^^^ 

Below we explicitly construct the MOOS for some particu- 
lar Lie algebra to be protected: 



1. For suppressing longitudinal relaxation of a single- 
qubit system along the z-direction, the MOOS contains 
only one operator cr,. 

2. For a general single-qubit system, two anti-commuting 
Pauli matrices, e.g., {cr,, cr,,), form an MOOS to protect 
all system operators H^l- 

3. For an L-qubit pure dephasing model in which H con- 
tains only the Pauli matrices along the z-direction, a 
choice of the MOOS is 

Note that only single-qubit operators are used. It is ob- 
vious that all the Pauli matrices in the set commute with 
each other There are totally L operators in the MOOS, 
which can be shown to be the minimum possible num- 
ber. 

4. For a general L-qubit system, to protect all operators 
of the system, i.e., the algebra su(2^), a choice of the 
MOOS is {cr?.",crf It is obvious that all the Pauli 
matrices either commute or anti-commute with each 
other Note that only single-qubit operators are used. 
There are totally 2L operators in the MOOS, which can 
be shown to be the minimum possible number. 

5. For a general M-level system, let us consider protection 
of all system operators which are diagonal in a certain 
orthonormal basis {\m}}^^i^. Suppose 2^"' < M < 2^. 
We can denote an integer number Q < m < M using a 
binary code as m = (m^ ■ ■ ■ m2m\) with m/ = or L We 
define a unitary Hermitian operator 

= 7-2 ^ |m><m|. (5) 

A diagonal operator of the form \m){m\ can be written 
in terms of {2;'') as 

T^ / + (-l)""Sf 
\m){m\ =11 . 

l=\ ^ 

Thus any diagonal operators can be constructed using 
{S!.''). Actually, by adding (2^-M) ancillary basis states 
{\m)]^^lj, the above-defined operators can be viewed 
as the single-qubit Pauli matrices along the z-direction 
({crl'*)) of an L-qubit system projected to the M-level 
subsystem. It is obvious that all operators in {E*'*) com- 
mute with each other Thus an MOOS is constructed. 
There are totally L such operators in the MOOS. Since 
an (L - l)-qubit system has at least (L - 1) operators in 
an MOOS and 2^-' < M < 2^, the MOOS {Sf con- 
tains the minimum number of operators (otherwise, we 
can construct a DD with less than M intervals to pro- 
tect the system to the first decouphng order, which is 
impossible according to Ref. 1I31I1 ). 
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6. For a given /, if M/2 is an integer, we can define unitary 
Hermitian operators anti-commuting with as 

4" = Yj {\m + 2'"' )<ot| + h.c.) , (6) 

m,=0 

which exchanges two basis states \m} and \m'} if m and 
to' differ at and only at the /th bit. Actually, the operator 
l};^ can be viewed as the Pauli matrix cr*'' of the fth qubit 
for a multi-qubit system. {li"|M mod 2' = 0) U {I.[\2' < 
M] forms an MOOS of the M-level system. 

Note 1 : The choice of an MOOS for a certain system is not 
unique. For example, in the MOOS for generating su(2^) of 
a general L-qubit system, the two Pauli matrices crl'* and cr*'' 
can be replaced with any two anti-commuting Pauli matrices 
of the fth qubit. 

Note 2: For a general M-level system, there may exist no 
MOOS to generate the whole su(M) algebra. For example, 
for an M-level system with M being an odd number, all oper- 
ators in an MOOS must mutually commute and therefore are 
all diagonal in a common basis. This MOOS of course cannot 
generate the whole su(M). Actually, if two operators Q and 
Q.' in an MOOS anti-commute, we have Tr(Q + Q'QQ') = 
2Tr(Q) - 0, which is impossible for odd M. 

For a general M-level system, the explicit operations given 
in Eqs. Q and (|6]l are in general difficult to implement in ex- 
periments. It is important to find a suitable set of operations 
for a given multi-level system. It should also be noted that 
above we have assumed that all the operators in an MOOS 
are protected to the same decoupling order In the CDD and 
NUDD schemes we will discuss later, different operators ac- 
tually can be protected to different orders. In those cases, the 
precision of protection of the whole Lie algebra is determined 
by the lowest decoupling order. 

B. Lowest order protection of system operators 

A general DD framework to protect a set of operators to the 
lowest order is the symmetrization procedure over an appro- 
priate DD group IHIqI]. Here we systematically give an ex- 
plicit scheme to protect system operators forming an MOOS, 
which facilitates the construction of higher order DD in the 
later part of this paper. 

Let us consider first protection of a single unitary Hermitian 
operator n. The Hamiltonian can be separated into two parts, 

H = Cn+An, (7) 

with 

Cn = {H + QHn)/2, (8a) 

An = {H-QHQ.)/2. (8b) 

Q. commutes with Cn and anti-commutes with An, i.e., 

OCnO = Cn, (9a) 
QAnO = -An. (9b) 



With an instantaneous control pulse Q applied at the middle 
of the evolution time, the evolution operator becomes 

UniT) = e-'"^'^ne-'"^'^ = ^^-inHnr/i^-.HT/i 

which commutes with Q. up to an error of 0{T-). In general, 
one may apply an additional pulse at the end of evolution so 
that 

[n]Uo(T/2)nUo{T/2) = (11) 

where Uq{t) = e"'^^ is the free evolution operator over time 
T, and the brackets around the operation at the end of the se- 
quence ([Q]) mean that the operation is optional. The effective 
Hamiltonian //n commutes with O up to 0{T). Therefore Q 
is protected to the first order. 

Following the method given above, we can preserve more 
operators {Q.^} in an MOOS by concatenation. The first level 
of control is 

Um = [Qi]f/o(r/2)Qif/o(r/2) = e-'"^^, (12) 

where the effective Hamiltonian H\ commutes with Qi up to 
0{T). By treating the effective Hamiltonian as a new Hamil- 
tonian on the second level, the propagator reads 

U2{T) = [Q2]f/i(7^/2)Q2t/i(r/2) = e-'"'-'', (13) 

where the effective Hamiltonian H2 commutes with Q2 up to 
an error of 0{T). And according to Theorem[T] H2 also com- 
mutes with Qi up to 0{T). A general first order scheme is 
achieved by using Eq. (fTTT i iteratively, 

Ul{T) = [^l]Ul-x{TI2)^lUl-i{TI2) = e-'"''\ (14) 

where the effective Hamiltonian Hi commutes with all oper- 
ators in the MOOS {Q/)f^j up to 0{T). Note that when the 
optional pulses are not used, only one pulse is applied at each 
time of operation. 

It should be pointed out that the current DD scheme con- 
structed by iteration is not the same as the symmetrization 
procedure described in Refs. fs)] and[l9l]. For M-level systems, 
when log2 M is not an integer, our scheme in general cannot 
achieve protection of all system operators. For example, when 
M is odd, our scheme can only protect operators diagonal in 
a certain basis. When M = 2^, all system operators can be 
protected by protecting the MOOS of size 2L. Several advan- 
tages of the MOOS-based DD, however, are worth mention- 
ing. First, there are systematic ways to construct higher order 
DD for protecting all operators in an MOOS and hence the 
Lie algebra generated from the MOOS (Sees. Ill CI and III Dl i. 
Second, for multi-qubit systems, our scheme automatically 
gives DD sequences involving only single-qubit operations 
(Sec. inil i. Third, there is an explicit scheme to incorporate 
pulses of finite amplitude to DD (Sec. lIVI i. Fourth, since the 
DD sequences on the inner levels use more control pulses than 
the sequences on the outer levels, qubits subjected to faster 
error sources can be protected on inner levels for economic 
use of control resources. For example, in a coupled electron- 
nuclear spin system, the electron spin, which has much faster 
decoherence than a nuclear spin does, should be controlled on 
the inner level. 
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C. Higher order protection by CDD 

To go beyond the lowest order protection, let us consider 
a general sequence of unitary operations (cr^) on the system, 
where cr^ can be the identity operator. The evolution from 
f = to r reads 



= 7"exp 



Jo 



{t)dt 



(15) 



where T is the time ordering operator, and //o-(f) = crj^ Hcri^ . 
for t € [tj, tj+i]. In the standard time-dependent perturbation 
theory formalism, the propagator is expanded up to the second 
order as 



:(T) = 1 + ^ /,,. + ^ hihj + ^ hij + 0(T\ 



(16) 



where hj = -/(fy+i - tj)Ha-(t) (with fo s and f„+i s T). 
In this paper we assume that this expansion converges, which 
means that the bath is bounded in spectrum. 

It was shown that if Uc{T) realizes the first order DD, i.e., 
[ZilLo'^i' Qj] - for ^ set of operators {Qj}, then the second 
order DD can be realized by the symmetrized evolution lUlIlt] 



Usdd{2T) = UcUc 

= l+2j]/, + - 2^/.,. 



I n \ 



(=0 



+ 0(T^), (17) 



V 1=0 / 



where Uq = e °e ' ■ ■ ■ e "-'e " is mirror-symmetric with Uc- 

For DD to even higher orders, in principle we can obtain 
the optimal sequences by solving Eq. ( fTSI l so that up to some 
order in the expansion, Uc{T) commutes with a given set of 
operators. The third order terms read 



i>j>k 



and in general i;j>;,>...>y„ //'//^ ■ ■■//'; n"=i ^ contains 
terms of the order Finding solutions becomes 

formidable when the DD order is high. 

If we are not concerned with the exponentially increas- 
ing number of control pulses, we can follow the idea of 
CDD lfl2ilT4ll to construct DD sequences in a systematic way 
to protect the operators to an arbitrary DD order The 

first order DD given by Eq. ( fT4l i is 2^ evolution operators 
U{T/2^) embedded in a sequence of control pulses {Q/). We 
denote this structure as 



U\^\T) =D.lUl-i{TI2)D.lUl-i{TI2) 



=Cn {U{T/2h} = e- 



iH\"T 



(18) 



which is the first order CDD as defined in Eq. ( fT4l i. The re- 
sultant first-order effective Hamiltonian commutes with 
the operators up to an error of 0{T). Here the se- 

quence Cn {■ ■ ■ ) makes the effective Hamiltonian commute 



with to a higher order. The evolution under the sec- 

ond order CDD 



t/f(r) = Cn{t/["(r/2^)} 



(19) 



is obtained by replacing the free evolution U{T/2^) in Eq. ( fTSb 
with U^j^\t/2^), an evolution operator under the first order 
CDD control. The resultant effective Hamiltonian com- 
mutes with up to an error of 0(T^). Iteratively, the Mh 
order CDD reads 



(T) = Cn 



[N-l] 



{T/2')] 



(20) 



which preserves any operators in up to 0{T^^^). Note 

that in the above construction of CDD, Theorem [T] has not 
been invoked, except in construction of the innermost level as 
in Eq. ( fT4b . The only requirement is that all the intervals of the 
outer level are equal so that the resultant effective Hamiltonian 
from the inner level of control is time-independent. 

An alternative construction of CDD is given as follows. We 
first construct a CDD sequence to protect Q.i up to an error of 
O (t^' ^ ' ) by recursion 



rC[0] 



C[/v,-i]^y,/2)Qif/^['^'-i](r/2), 



(21) 



with U^^ (T) = U{T) and the superscript C denoting the nest- 
ing scheme of CDD. By defining U^^^^\T) = U^^'\t\ 
we can construct a further level of CDD to protect Q2 up to 



O {t'^-^^^, by the recursion 



^[^■;^-'l(r/2)Q2t/;^'^'f-"(r/2). (22) 



Similarly, we have the propagator by recursion 



jjC{NuN2,...,N,fi] ^ ttCIN,,N2,...,Ni].j,. 

TjC[Ni,...,NL](r^. _n jtC[Ni,...,Nl-i,Nl-1] 



'ni,n2,...,nt 



(r/2)Qi 



TTC[Ni,...,NL-uNL-l](r /l^ 



(23a) 



(23b) 



According to Theorem [T] the inner levels of DD are un- 
affected by the outer levels of control. Thus the evolution 
(T) commutes with Q.\,Q.2, . . . ,Q.l^o the orders 



C[Nx,N2,...,Nl] 

ni,n2,...,nz. 



U\ 

N\,N2, ■ ■ ■,Nl, in turn. An advantage of this construction is 
that the errors induced by Hi are eliminated independently 
and to different orders {Ni}, which allows protecting operators 
with stronger error sources to higher orders. For example, 
usually for spin qubits under strong external magnetic field, 
the pure dephasing is much faster than the population relax- 
ation, so it is favorable to protect the phase correlation on the 
inner level and to a higher CDD order. 

The number of operations required in preserving the oper- 
ators is ~ 2'^^ for the CDD scheme in Eq. ^ or ~ 2^^' 
for that in Eq. ( |23] ). They increase exponentially with the DD 
order. Even though the exponentially increasing number of 
control pulses does yield significant improvement of precision 
(through reduction of the coefficient in front of the power of 
time r'^+i) iUllllll], implementation of CDD to high or- 
ders is challenging in experiments since errors are inevitably 
introduced in each control pulse. 
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D. Higher order protection by NUDD 

If there is only one unitary Hermitian operator Q. to be pre- 
served, operations of Q applied at the UDD timing lfl6ll 



T„ = T sin^ 



2N + 2 



, for « = 1 , . . . , A^, 



(24) 



during the evolution from Tq s to T^+i = T protect the 
physical quantity D. to the A^th order IjtIi . Explicitly, the prop- 
agator under control 



U^^^^T) = Q^C/o(r^) ■ • • Qf/o(Ti)Qf/o(To), 



(25) 



commutes with f2 up to an error of 0(T^^^). Here the evolu- 
tion intervals are 



T„ = T„ 



T 



nn 



cos ■ 



■ cos ■ 



(« -I- l);r 



1 



(26) 



When there are more than one unitary Hermitian operators 
{Q/)^ I to be protected, the question is whether we can con- 
struct NUDD so that the number of control pulses scales 
polynomially with the protection order. A known example 
is QDD, in which UDD of, e.g., cr, and cr, are nested. For 
general cases, we establish the following theorem as the basis 
of NUDD. 



1. A theorem on UDD control of time-dependent systems 

Theorem 2. For a finite-norm time-dependent Hamiltonian 
H(t) defined in [0, T\ an Nth order UDD control with N op- 
erations of unitary Hermitian operator O applied at Ti, T2, 
■ . ., T]\f preserves Q up to an error o/(9 ^T^^'), ;/ 



H(T„ + ST„) = H(T„^i - ST„) = HisTi), 



(27) 



for s e [0, 1] and t„ — T„+i — T„, i.e., the Hamiltonian has the 
same form as a function of the relative time between adjacent 
operations and is symmetric within each interval. 

Note: The previous extension of UDD to time-dependent 
systems requires that the Hamiltonians be analytic (having 
smooth time-dependence) IztI l29i [stIi . In Theorem |2] the 
Hamiltonians are not required to be analytic but with cer- 
tain symmetries. The symmetry requirements on the time- 
dependence of the Hamiltonians can actually be fulfilled by 
designing the timing of DD sequences on the inner levels so 
that recursive nesting of DD is possible. 



Proof. The evolution under the control of reads 
U(T) = Q.'^VN^yN-i ■ ■ ■ QViOVo, 



(28) 



with the evolution operator 
y„ sTexp 

= Te exp 



XT„+\ 
H{t)dt 



\6)d9 



(29) 



where Te stands for ordering in 6, and 

A^-i- 1 



H''\e) 



-H(t), 



(30) 



where = ^ + for t e (T„,T„+i]. The symmetry 

requirements given in Eq. ( |27] ) are transformed to 



0) = H"'\^^^^7T -6)^ W\e). (31) 



N+l 

The Hamiltonian H'^^{0) can be separated into two parts. 



H'^\e) = Cie) + A{e), 

C(0) = [//''='(f) + Qi/"='(0)Q] /2, 

A(0) = [//'^'(e) - nH'''\e)n\ 12. 



with 



(32) 

(33a) 
(33b) 



C(0) and A{ff) commute and anti-commute with the operator 
Q., respectively. 

Now we rewrite the propagator as 



f/(r)=Teexp 
where 

G(0) = 



-iT ^ G{ff)[C(0) + F{e)A{ff))d0 



1 



(« + l)7t 



N-h 1 



F(0) = (-!)«, 



(34) 

(35a) 
(35b) 



for e {^,^-^1 Thus, the part of Hamiltonian C{0) 
that commutes with O is modulated by the step function G{0) 
which has step heights given by the UDD intervals, and the 
part of Hamiltonian A{0) that anti-commutes with Q. is modu- 
lated by G{0) and the periodic modulation function F{0). Fur- 
thermore, both C{0) and A{0) have the same symmetries as 
//■^•^'(fl) in Eq. dSTT i. The symmetries of the time-dependent 
Hamiltonian and the modulation functions make them have 
particular Fourier expansions, which lead us to a proof of 
the theorem in a procedure similar to the proof of UDD in 
Ref. iH. 

The Fourier expansions of the modulation functions and the 
time-dependent Hamiltonians are 



F(0) 



G(0) = 2§isin[2^(A?+l)0+0], 

k=0 

CO 

^/,sin[(2fc+l)(A^+l)0], 

A-=0 

CO 

Y^ck cos[2kiN -¥ 1)0], 

k=0 

CO 

A(0) = ^ at cos[2;t(A^ H- 1)0]. 



C{0)^ 



(36a) 
(36b) 
(36c) 
(36d) 



k=0 



Here the operators and commute and anti-commute with 
Q, respectively. The features of these Fourier expansions to 
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be used in the proof below are: (i) Both C and A contain only 
cosine harmonics of order of even multiple of (A^ + 1); (ii) 
F{9) contains only sine harmonics of order of odd multiple 
of (A^ + 1); G{0) contains only sine harmonics of an order 
differing from an even multiple of (A^ + 1 ) by + 1 or - 1 . 
With the product-to-sum trigonometric formulae, we have 



t/(r) = reexp 



-iT^ {c{e) + A{e))de 



(37) 



with 



c{e) =G{e)C(e) = ^ q fim[2k(N + \)e ± e], (38a) 



A{e) =G{6)F(0)A{0) 

= ^ a* cos[(2fc + 1){N + 1)0 + 0]. 



(38b) 



A straightforward method is to expand U{T) according to 
the standard time-dependent perturbation theory. It should be 
noted that such perturbation-theoretic expansion requires that 
the modulated Hamiltonian have bounded norm. In the expan- 
sion, the terms which do not commute with Q must contain an 
odd times of {5^) (since Q. anti-commutes with {5,^)). The ex- 
pansion coefficients can be written as 



{-iT)"x 

Jo Jo " Jo 



)d0i---d0„, 
(39) 



with yl'^{0) = sin[2k(N + 1)0 ±0] associated with an operator 
Q, and yl'^(0) = cos[(2A; + l)(N + 1)0 + 0] associated with 
an operator at, for aj 6 {c, s) and rjj € {+, -}. By induction 
and repeatedly using the product-to-sum trigonometric formu- 
lae, one can straightforwardly verify that the coefficients in 
Eq. (39[ vanish for n < N and y^'^ appearing an odd number 
of times. Thus vanish any terms in the expansion which con- 
tain products of an odd number of operators in [cik] and have 
a power of T lower than (A^ -i- 1). 

□ 



2. NUDD 

For a time-independent Hamiltonian H under DD control 
of instantaneous operations of a unitary Hermitian operator 
Q. applied at fi, f2, . . . , tt^t, the evolution U{t) from fo = to 
fiV'+i = T is equivalent to the evolution under a time-dependent 
Hamiltonian 



(40) 



for f 6 (f„,f„+i]. Such time-dependence is not analytic. If 
A^' is an even number and the operation sequence is symmet- 
ric, the time-dependent Hamiltonian H{t) is time symmetric 
in [0, r]. Thus, according to Theorem [2] in a UDD sequence 
of an operator Q applied at T^To, ■ ■ between To = 



and Tn+i = T, each interval of free evolution e-'fiiT„+i-T„) ^^j^ 
be substituted with the evolution inserted by a sequence of 
another operation Q' applied at r„.i, T„2, ■ ■ ■ , T„^n' between 
T„fi = T„ and T„^n'+\) = T„+\, with the same symmetric struc- 
ture in all intervals, i.e.. 



Tn+l - Tn,N'-k 


- T„,k+i 


-T„, 


(41a) 


Tn,k - T„,k' 




Tm,k' 


(41b) 


T,i+i - T„ 


Tm+ 1 


- T ' 



In particular, the inner level control of Q.' can be chosen as an 
even order UDD. 

Now we describe the construction of NUDD for protecting 
a set of unitary Hermitian operators First, the A^^th 

order UDD sequence of is constructed with pulses applied 
at 



T„, = T sin^ — '-^^ — , 
2Nl + 2 



(42) 



between - and T^^+i = T as the outermost level of con- 
trol. Nl could be either odd or even. Then the free evolution in 
each interval is substituted by the Nl- i th order UDD sequence 
of Hl-i applied at 



T - T 



-{Tnt+i - T'nz.) sin- 



2NL-1+2' 



in each interval between T, 



nt,0 



T„^ and T„^^Nl- 



(43) 



with Ni-\ being an even number So on and so forth, the 
/th level of control is constructed by applying Ni times of 



Q/ in each interval between T, 



^ni,...,nM,N,+\ 



iiL....,niti,0 



= T, 



;iL,...,n/+2,/i;+ 



and 



= T, 



Hz,,...,n/+2,n/+i + 



1 at 



T - T 



+ [r„,,...,„„,,„,„+i - r,„,...,„,^,] sin^ (44) 

with A^; being an even number We denote the evolution un- 
der such NUDD as U^^^^i'^--'I^''\T), where the superscript U 
denotes the nesting of UDD sequences. 

According to Theorem |2l the outer levels of UDD control 
are not affected by the inner levels of even-order UDD control. 
And according to Theorem [1] the inner levels of DD control 
are not affect by the outer levels of control since {Q;) is an 
MOOS. Thus each operator Hi is protected up to an error of 
O (r'^'^^y The number of control intervals is 



N' 



V[Ni .N2-;Nl] 
pulse 



^(Ni + 1){N2 + 1)---{Nl + 1), 



(45) 



increasing polynomially with the decoupling order. 

III. NUDD OF MULTI-QUBIT SYSTEMS 

A. General multi-qubit systems 

To protect a multi-qubit system to a given order of preci 
sion, we just need to protect the operators in an MOOS de 
scribed in Sec. Ill Al by using the method depicted in Sec 
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for the first order preservation, or by using the NUDD scheme 
in Sec. IIIDI for higher order preservation. Since for general 
muhi-qubit systems, the MOOS can be chosen as a set of 
single-qubit operators such as the Pauli matrices {cr*.'\ cr''')^j, 
NUDD can be implemented with only single-qubit flips. 

In suppressing the relaxation or pure dephasing to the 
first order, i.e., in protecting the MOOS {erf or {o-f 
the first order scheme in Sec. IIIBI requires 2^ pulse inter- 
vals. For suppressing decoherence in general cases, the 
first order scheme requires 4^ intervals to protect the MOOS 
Such numbers of intervals are actually the min- 
ima required for protecting L-qubit systems, as proven in 
Ref. Bill . In preserving the coherence of a multi-qubit system 
to an arbitrarily high order of precision by NUDD, the num- 
ber of intervals given in Eq. ( |45] ) increases polynomially with 
the decoupling orders, much less than that required in CDD. 
A question is whether NUDD is optimal or nearly optimal in 
terms of the number of control pulses. For QDD control of 
one qubit, numerical check up to the fourth order indicates 
that NUDD is nearly optimal, differing from the optimal so- 
lutions by less than 3 control pulses 112611 . NUDD of a larger 
MOOS, however, can be shown to be far from the optimal 
in the second decoupling order; For an L-qubit system suf- 
fering pure dephasing, an L-level NUDD in the second order 
requires ~ 3^ control intervals, while the SDD, which uses 
two mutually symmetric first order DD sequences to realize 
the second order control (see Sec. IIICI ), requires only 2x2^ 
intervals. We expect that in higher orders of DD, there exist 
DD schemes (using only single-qubit control for multi-qubit 
systems) much more efficient than NUDD. But no explicit so- 
lutions are known to us, except for a few numerical solutions. 

B. Discussions on QDD 

In particular, for a single-qubit system, NUDD reduces to 
QDD with even order DD on the inner level. Explicitly the 
nested sequence U]^^^^^''^-\T) is the QDD sequence protect- 
ing the MOOS operators and a-^ to orders 2N\ and A^2, 
respectively. Thus based on Theorems[T]and |2] the validity of 
QDD with even order UDD on the inner level is proven. 

But Theorem|2]does not apply to the case of odd order UDD 
control on the inner level. Actually, when the inner level UDD 
has an odd order, which breaks the symmetry condition of the 
theorem, the outer level UDD may be spoiled. For a specific 
example, let us consider the control of a Hamiltonian like 

where Jq, Ji, J2, and 7i_2 are arbitrary bath operators, and Qi 
and Q2 are two system operators forming an MOOS (such as 
Qi = and Q.2 = for a single-qubit system). We choose 
the inner level control as the first order UDD of Oi and the 
outer level as the second order UDD of Q.2. The propagator of 
this NUDD is 

x(pie-'"'Q.ie-'"^), (46) 



where the evolution in each pair of parentheses corresponds 
to the UDD control on the inner level and r = r/8. The time 
expansion gives 

f/J^['n^'(r) = C/b + QiO{T^) + QiQ2C»(r2) ^ q^j^3^^ ^47^ 

where [/^ is a pure bath evolution operator. Thus even though 
a second order UDD sequence of Q.2 is applied, Q2 is pre- 
served only to the first order. 

The above example indicates that the effective Hamiltonian 
resulting from the inner level UDD control can not be writ- 
ten into an analytic form. Otherwise, according to Ref. IztIi . 
which establishes the performance of UDD on analytically 
time-dependent systems, the outer level UDD should not be 
affected. As proposed in Ref. IztIi . for an NUDD evolution 

U"^%f\T) = Q?'f/^|^'l(r^)Q2 ■ ■ ■ <^'l(rOQ2f/r^(ro), 

one can define the effective Hamiltonian resulting from the 
inner level of control as 

//eff(r„ + r) = / [d.U^'^'^t)] [<^'1(0]"'" , (48) 

for t € (0, T,,]. The outer level UDD can be viewed as act- 
ing on this effective Hamiltonian. The Hamiltonian defined 
in Eq. (|48j, however, is only piecewise analytic and is even 
discontinuous at T„'s. Therefore, the theorem established in 
Ref. ll27ll about UDD control of analytically time-dependent 
systems does not apply to QDD. 

Thus, the complete proof (or disproof) of the validation of 
QDD is still an open question. We should mention that if QDD 
is proven valid for one qubit, according to our formalism of 
nested DD, the same NUDD is also valid for any two opera- 
tors forming an MOOS. In this way the QDD control can be 
generalized, in particular, to the protection of two-qubit sys- 
tems from dephasing and disentanglement, etc. 

IV. DD BY PULSES OF FINITE AMPLITUDE 

The ideal instantaneous pulses are not realistic in experi- 
ments since they would contain an infinite amount of energy. 
Of course, when the pulses are mush shorter than the other 
timescales of the system and the bath, it is a good approxi- 
mation to treat them as infinitely short. But if this condition 
is not satisfied, it is of interest to consider DD by pulses of 
finite amplitude. The problem of first order DD with finite- 
amplitude pulses has been considered within the Eulerian DD 
framework |[38ll and in a geometric picture ll39ll . It is possible 
to achieve arbitrary control precision by recursive construc- 
tion of pulse shapes |40]. DD of single-qubit systems using 
finite-amplitude pulses up to a control error in the second or- 
der of pulse durations has been presented in Refs. ll32i liUl . 
In UDD, finite-amplitude pulses of higher orders of control 
precision can also be incorporated ll20ll . 

Here we consider the general case of DD by finite- 
amplitude pulses. Let us consider a short-pulse operation by 
the Hamiltonian 

Hnit) = v{t)Q, (49) 
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where Q is a unitary Hermitian operator. We aim to design 
the pulse shape of v(f) such that the evolution during the pulse 
control approximates the ideal 5-pulse control up to a certain 
order of the pulse duration Tp, i.e., 



Jo 



[H + Hnit)] dt 



= e-"^''^-''''"Pae-"'" + 0{Tp" ), (50) 

where Pq = exp |^-/ j^'' Hn{t)dt^ is the desired instantaneous 
control applied at the time r,. In particular, we need the 
n pulse Pn = Q up to a trivial global phase factor Un- 
fortunately, a no-go theorem established in Ref. ll32l l4lll re- 
stricts that instantaneous n pulses can not be approximated by 
a finite-amplitude pulse with error lower than O (t^) without 
perturbing the bath evolution. Thus, here we focus on the first 
order pulse shaping with Mp = 2. 
We write the evolution operator as 



where 



Un -T exp 



Jo 



-I I Ha{t)dt 
Jo 



(51) 



(52) 



with Hnit) = e''^^''''''>Hn{f)e'''^^'''''\ The correction term is 
hnit) = Hnit) - Hnit) 
^ it-T,)" 

k=\ 



vit)J]- 



k\ 



[iH,[iH,---[iH,D] ■■■]]. (53) 



k folds 



(54) 



We want to design y(f) such that the control error 

6Pa =Un-Pn 

The leading order term in 6Pn is 

= T l^e-'fo" j\t - Ts)vit)[H, Q]^fr| 

it - T,)v(r)e-''.t " ""('1d'[H, Q]e-'/o ""<''1'''dt. (55) 



-f 

Jo 



Using Q - 1 and Hnit) - vit)D., we have 



f " Hnis)ds 

e •>•! - cos 



r 

Ji\ 



vis)ds 



;0 sin 



I \'is)ds 
Jt, 



(56) 



Now we decompose the Hamiltonian H into two parts as H - 
A + C, with A and C anti-commuting and commuting with Q, 
respectively. The leading order error term in Eq. dSSl l becomes 



/7^" = [A,n]77„-/Q[A,Q]?7i2, 



(57) 



where 



711= r it -T,Mt) cos [4>o-i^it)]dt, (58a) 
Jo 

'712 = r it- T,)v(f) sin [4>Q - iA(0] dt, (58b) 
Jo 



with tffit) = '2. f' vis)ds and = JJ" vis)ds - J^' vis)ds. To 
eliminate the leading order error, we just need to make ?7i i = 
?7i2 = 0, which are the same as those derived in Ref. Il32ll for 
single-qubit flip control. 

Using the finite-amplitude pulses designed as depicted 
above, we can realize DD based on an MOOS up to an error in 
the second order of the pulse duration. Note that in some DD 
schemes, such as CDD U^\'^''-^f''\T) as shown in Eq. ( |23] |, 
operations of different Q/'s may coincide. For example, Qi, 
Q.2, and Q3 coincide at the end of the sequence U^'^^^^^^^\T). 
In this case, we can define a new unitary Hermitian operator 
Q' as the product of the operations (such as QiQa^a in the 
example), and design the pulse Hn'it) - vit)Q.' to achieve 
the operation up to an error of (9(t^,). Such finite-amplitude 
pulse operation, however, would involve multi-qubit interac- 
tions and may not be easy to be implemented in experiments, 
unless Q.' happens to be a single-qubit operation. This prob- 
lem, fortunately, does not exist in NUDD, since there no two 
operations coincide, which stands for another advantage of 
NUDD over CDD. 

Ref. 1I20I1 has presented a method to implement UDD by 
higher order shaped pulses. At the first sight, it seems that 
those pulses can be incorporated in NUDD. However, the 
method in Ref. 1I20I1 requires a starting and a stopping pulse 
to protect the UDD sequence. In NUDD, such starting and 
stopping operations will be mixed up with the outer level con- 
trol applied at the same time. Then the operations on the 
outer levels need to be redesigned, which may be much more 
complicated than the design in Ref. ll20ll since different op- 
erations may interfere with each other and multi-qubit inter- 
actions may be involved. Explicit implementation of NUDD 
with finite-amplitude pulses of higher order control accuracy 
is an interesting topic for future work. 



V. CONCLUSIONS AND DISCUSSIONS 

Based on two theorems, we have presented explicit 
schemes of dynamical decoupling to preserve operators in an 
MOOS (i.e., unitary Hermitian operators which either com- 
mute or anti-commute with each other) to arbitrary decoupling 
orders for quantum systems arbitrarily coupled to quantum 
baths. All system operators in a Lie algebra generated from 
the MOOS by commutation, anti-commutation, linear combi- 
nations, and repetitions are also preserved. Theorem [1] states 
that the inner levels of DD control are unaffected by the outer 
levels if the control operations are elements of an MOOS. 
Theorem|2]states that UDD still works if the Hamiltonians in 
different intervals have the same function form of the relative 
time and are symmetric, regardless of the analytic properties 
of the Hamiltonians. These theorems enable a construction of 
higher order DD by nesting UDD sequences of even orders. 
NUDD protects system operators in a Lie algebra generated 
from an MOOS to an arbitrary order of precision. For multi- 
qubit systems, any physical quantities can be protected, and 
NUDD can be implemented by single-qubit operations. For 
single-qubit systems, NUDD reduces to QDD with even order 
UDD on the inner level. Thus the theorems provide a rigorous 
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proof of the validity of QDD with even order DD on the inner 
level. 

NUDD achieves a desired decoupling order with only a 
polynomial increase in the number of pulses, with exponen- 
tial saving of the number of pulses as compared with CDD of 
the same decoupUng order. In suppressing the general deco- 
herence, the number of pulses still scales exponentially with 
the number of qubits. Such exponential increase, indeed, is 
required by a theorem which sets the minimum number of 
control intervals to be 4^ or 2^ for protecting a general or 
pure dephasing L-qubit system to the first decoupling order, 
respectively IBlll . 

For Hamiltonians of certain structures, such as the Hamil- 
tonians of qudit systems with bipartite interactions, reduction 
in the number of pulses is possible ll42ll . Fewer levels of 
nesting are required if the structures of the Hamiltonians are 
exploited and a proper MOOS is designed. The number of 
pulses can also be greatly reduced if we protect only some 
logically encoded qubits or some particular states of the sys- 
tem IHIIllllillll. For example, if we choose the MOOS as 
{cr'/* igi c/f* , (T^P, crp'} for a two-qubit system, the only possi- 
ble noise generator Issll after protection is cr*'' ® cr?*, which 
commutes with all the elements in the MOOS; in this case, the 
logical qubit a\ TT) + P\ ii) or a'\ tJ.) + 13'\ XT) is protected 



by only three levels of nesting. 

NUDD protecting two operators forming an MOOS is near- 
optimal and has the same timing as QDD. In general, how- 
ever, NUDD is by far not optimal, since a large nesting level 
L requires much more control intervals than the symmetrized 
DD for achieving the second decoupling order. An interesting 
question for future study is how to construct optimal or nearly 
optimal higher-order DD for general multi-qubit or multi-level 
systems. 

For realistic implementation of DD, we have derived the 
conditions for finite-amplitude pulses to simulate ideal opera- 
tions up to an error in the second order of pulse duration, and 
the conditions reach the same results as for single-qubit flip 
control given in Ref. I.3Z1 . Thus we can apply the pulses de- 
signed in Ref. lf32ll to the higher order DD schemes for general 
quantum systems. 
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